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Abstract 

We present the first calculations of a symmetry conserving configuration mixing method (SCCM) using time-reversal 
symmetry breaking Hartree-Fock-Bogoliubov (HFB) states with the Gogny D1S interaction. The method includes 
particle number and tridimensional angular momentum symmetry restorations as well as configuration mixing within 
the generator coordinate method (GCM) framework. The nucleus 32 Mg is chosen to show the performance and reliability 
of the calculations. Additionally, 0^, 2+ and 4+ states are computed for the magnesium isotopic chain, where a noticeable 
compression of the spectrum is obtained by including cranked states, leading to a very good agreement with the known 
experimental data. 
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A trustworthy description of the spectra of the atomic 
nuclei is one of the main goals of low-energy nuclear struc¬ 
ture theory. The interacting shell model [1, 2] is likely 
the most widely used and successful tool to compute ac¬ 
curately low-lying spectroscopic properties. However, shell 
model applications are limited to regions not far away from 
shell closures where manageable valence spaces can be de¬ 
fined. On the other hand, microscopic self-consistent mean 
field methods (SCMF) [3] based on nuclear energy den¬ 
sity functionals such as Skyrme, Gogny and/or Relativistic 
Mean Field (RMF) can be in principle used throughout the 
whole nuclear chart. In order to apply these methods to 
the study of nuclear spectra, they have to be extended by 
including beyond-mean-field (BMF) correlations. In par¬ 
ticular, symmetry conserving configuration mixing meth¬ 
ods (SCCM) are the most natural extensions of SCMF 
approaches and have shown a fair performance in describ¬ 
ing qualitatively not only nuclear spectra, but also ground 
state properties, electromagnetic transitions and decays. 
Unfortunately, quantitative accurate predictions have not 
been reached so far, mainly due to the lack of time-reversal 
symmetry breaking intrinsic states within the existing im¬ 
plementations of the SCCM methods. In this letter we 
present an extension of the SCCM framework, based on the 
generator coordinate method (GCM) with particle num¬ 
ber and triaxial angular momentum projections, that in¬ 
cludes cranking intrinsic states. In the numerical appli¬ 
cations we use the finite range density dependent Gogny 
interaction (D1S parametrization [4]). In the earliest im¬ 
plementations of the SCCM method only axially symmet- 
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ric intrinsic states were considered [5-9]. A major break¬ 
through towards a better description of the nuclear spectra 
within the SCCM framework was the inclusion of the tri¬ 
axial degree of freedom [10-12]. A further step forward 
was the first implementation of a SCCM method based on 
a Skyrme pseudo-potential for odd nuclei by B. Bally who 
obtained very promising results in describing the bench¬ 
mark nucleus 25 Mg [13]. 

On the other hand, the angular momentum projection 
with the energy density functionals mentioned above is 
performed after the energy variation. Therefore, the con¬ 
sideration of only intrinsic wave functions with zero an¬ 
gular momentum content (< J x >=< J y >=< J z >= 0) 
in the the current SCCM calculations tends to favor the 
ground states with respect to other excited states and a 
stretching in the spectra is usually found with respect to 
the experimental values. The addition of time-reversal 
symmetry breaking intrinsic states (< J x >7^ 0) obtained 
by the cranking procedure will thus increase the varia¬ 
tional space for excited states and will provide a better 
description of the spectrum. Pioneering angular momen¬ 
tum projection of cranking states have been reported with 
schematic pairing plus quadrupole interactions [14-16] and 
with Skyrme energy density functionals [17, 18]. However, 
neither configuration -shape- mixing nor, in the case of 
Sky me interactions, pairing correlations, were taken into 
account. The aim of this letter is to present the first 
results of the extension of the SCCM method described 
in detail in Ref. [11] (and references therein), including 
now time-reversal-symmetry breaking intrinsic states in¬ 
troduced through cranking calculations. 

The starting point is the construction of a set of intrin¬ 
sic many-body states having different deformations and 
intrinsic angular momentum. Such states, |/?2, 7 , J c ) = 
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Figure 1: (color online) (a)-(c) PN-VAP and (d)-(f) Particle num¬ 
ber and angular momentum projected potential energy surfaces for 
different values of the cranking angular momentum J c for the nu¬ 
cleus 32 Mg. Gogny D1S interaction is used here. The contour plots 
are separated in energy by 1.0 MeV. Each PES is normalized to the 
energy of their corresponding minima, i.e., (a) -249.902 MeV, (b) - 
247.910 MeV, (c) -246.789 MeV, (d) -252.924 MeV, (e) -252.021 MeV 
and (f) -250.463 MeV. The black bullets in (a)-(c) are the states in¬ 
cluded in the GCM calculation while the yellow squares in (c) are 
the states analyzed in detail in Figs. 2- 3. 


|), have the structure of Hartree-Fock-Bogoliubov (HFB) 
states [19] and are found by minimizing the particle num¬ 
ber projected HFB energy 1 , i.e.: 

E 'j c {P2,l) = ^JpNpzy^UJciJx) - Aq 20 (Q2o)-A g22 (<322) 

( 1 ) 

where P N ( Z ) is the neutron (proton) number projection 
operator [19]. This is the so-called variation after particle 


1 For the sake of simplicity, we will write down throughout the 
text any energy kernel as an expectation value of a hamiltonian op¬ 
erator. However, Gogny interactions contain a density-dependent 
term which prevents such a notation rigorously [20]. Nevertheless, 
this term is handled separately in a similar fashion as in Refs. [6, 11], 
and the following notation can be still used. 


number projection method (PN-VAP) [21]. The first term 
in the r.h.s. of Eq. 1 is the particle number projected en¬ 
ergy and the last terms correspond to the constraints on 
the cranking angular momentum J c and on the quadrupole 
deformation of the system (/? 2 , 7 )- The Lagrange multipli¬ 
ers co>j c , A q20 and X q22 ensure the conditions: 

(Jx) = VJc(J c + 1 ); (Q20) = #20; (Q22) = Q22 ( 2 ) 

where J x is the x-component of the angular momentum 
operator and Q^ix with fi = - 2 ,- 1 , ..,2 is the /i component 
of the quadrupole operator. The deformation parameters 
mentioned above are defined as: 

P2 cosy P2 sin 7 /~5~ 4tt 

q 20 -—^' q 22 ~^/W' C ~\^ (3) 

being A the mass number and tq = 1.2 fm. In the present 
work we have imposed the parity conservation as a self- 
consistent symmetry of the intrinsic states: V\) = |), be¬ 
ing V the parity operator. Therefore, neither negative 
parity states nor odd-multipole deformation -such as the 
octupole- degrees of freedom are explored here. Further¬ 
more, these states are invariant under the so-called simplex- 
x, S x \) = Ve~ l7V ^ x \) = |), and the T-simplex-y, Sy\) = 

V e ~ Z7rJy T\) = |) symmetries, being T the time-reversal 
operator. The last condition is chosen to have real co¬ 
efficients in the HFB transformation, and the simplex-x 
symmetry is very suitable to perform cranking calcula¬ 
tions (Eq. 1). The set of operators {V, S x , Sy} are the 
three generators of a subgroup of the more general point 
group Dj h [22, 23]. The latter has an additional gener¬ 
ator, e.g., the time-reversal operator. Although the use 
of self-consistent symmetries constrains the inclusion of 
correlations within the mean-field approach through the 
spontaneous symmetry-breaking mechanism, they are im¬ 
posed to reduce the computational burden. In the present 
case, we will also exploit these self-consistent symmetries 
to test the performance of the method since they provide 
non-trivial checks that help to identify possible inconsis¬ 
tencies. For instance, the choice of the collective coordi¬ 
nates (^ 2 , 7 ) divides the possible quadrupole deformations 
in six sextants, depending on the range of the angle 7 [19]. 
As a result, the values of 7 equal to 0°(60°), 120°(180°) 
and 240° (300°) correspond to prolate (oblate) axial defor¬ 
mation and they are related by the different orientations of 
the principal axes of inertia with respect to the z-axis [19]. 
If J c = 0, the intrinsic wave functions do not break the 
time-reversal symmetry and the energy is independent on 
the orientation of the axes, being all of the sextants com¬ 
pletely equivalent. However, if J c 7 ^ 0, the energy will 
depend on the orientation of the principal axes of iner¬ 
tia with respect to the intrinsic rotation axis, in our case, 
the x-axis (see Eq. 1). Therefore, the intrinsic states are 
only invariant under the subgroup of the D\ h point group 
mentioned above and the sextants are now symmetric only 
with respect to the 7 = (120°, 300°) direction. 
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We check this symmetry by performing PN-VAP cal¬ 
culations in the (/? 2 ,0° < 7 < 360°) plane for three values 
of the cranking angular momentum J c = 0 (time-reversal 
symmetry conserving), 2 and 4, selecting the nucleus 32 Mg 
as an example. In Fig. l(a)-(c) we show such potential en¬ 
ergy surfaces (PES) -first term in the r.h.s. of Eq. 1 . Here, 
the intrinsic states were expanded in nine major spherical 
harmonic oscillator shells and the number of points in¬ 
cluded in the mesh of each PES is 502. We notice first 
the equivalence between all of the sextants in the case 
where the time-reversal symmetry is preserved (J c = 0, 
Fig. 1(a)). Such a redundancy is reduced to the half of the 
plane separated by the 7 = (120°, 300°) axis for J c = 2 
and 4 (Figs. 1(b)- (c)) as expected. We find the absolute 
minimum of the J c = 0 PES in the spherical point as 
it is presumed from the neutron magic number N = 20. 
The energy grows more rapidly along the oblate directions 
than in the prolate ones. Additionally, a second minimum 
-around 1 MeV higher- is obtained at axial prolate configu¬ 
rations with /? 2 = 0.45. For larger values than /3 2 ~ 0.7 the 
energy increases quickly also along the prolate lines. For 
J c = 2,4 the PES resemble the J c = 0 one except for the 
values along 7 = 120 °, where the energy is not as favored 
as in 7 = 0° and 7 = 240°. The minima of these surfaces 
appear at such prolate configurations with /? 2 = 0.45, as in 
the J c = 0 case but shifted to higher energy values, around 
« 2 MeV and « 3 MeV for J c = 2 and 4 respectively. 

The intrinsic many-body states, |), break also the ro¬ 
tational invariance of the hamiltonian and these quantum 
numbers can be restored by projecting onto good number 
of particles (N,Z) and angular momentum (J, M): 


| JM; NZ; a ; ft, 7, ft) = £ Je) | JMK; NZ; ft, 7, ft) 

K 

(4) 

where K = — J, — J + 1,..., J is the component of the an¬ 
gular momentum in the body-fixed 2 -axis and the states 
given in the r.h.s. are defined as: 

\JMK-NZ-P 2 , 1 ,J c )=P j mk P n P z \P 2 , 1 ,J c ) (5) 

being P^k th e angular momentum projection operator [19]. 
Additionally, the coefficients 7 Jc ) and the projected 

energies Ep° Jc are found variationally by solving the 
so-called Hill-Wheeler-Griffin (HWG) equations in the K- 
subspace [19]: 



j-nz 

KK'(/3 2n ,J c ) 


zpJcr K[J\NZ 
' KK'(f3 2 ,j,J c 



= 0 


(6) 

where a = 1, 2,... labels the possible solutions of such gen¬ 
eralized eigenvalue problems and: 
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Figure 2: (color online) Probability distributions of projections K 
-VF^(/32,7, Jc)~ for even (left panel) and odd (right panel) values of 
the angular momentum J for the intrinsic states: (a)-(b) (^ 2 , 7 , Jc) = 
(0.5,10° A) an d ( c )-(d) (/#2»7 ,Jc) = (0.5, 230°, 4). Distribution of 
probabilities of (e) even values and (f) odd values of the angular 
momentum J -W J ((32 : r y, J c )~ for the same intrinsic HFB-type wave 
functions. 


are the norm and hamiltonian overlaps. From the above 
definitions, we can obtain some useful properties of the 
state P N P z |) such as [ 16 , 19 ]: a) the probability distribu¬ 
tion, ^(ft,7,</ c )i of finding an eigenstate of the angu¬ 
lar momentum | JK); b) the total probability distribution, 
W J (fa, 7, Jc) 1 of finding a value of the angular momentum 
J; and c) the projected energy ^^(^2,7, J c ). 


wycft,7,ft) = 


N j ’ nz 

Jy KK(p 2 ^,J c ) 


<ft, 7 , ftiP^lft^, ft) 


( 9 ) 


W 7 ft, 7 ,ft) = ^Wjhft, 7 ,ft) 

K 


K (ft,-7Vc) 
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J\N Z 

KK(p 2 , 7 ,J c ) 




J;NZ 

KK(p 2 ,-iJ c ) 


( 10 ) 

( 11 ) 
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Figure 3: (color online) Particle number and angular momentum 
projected energy overlaps E^(/32, 7 , J c ) as a function of K for the 
intrinsic wave functions (a) (/?2,7, Jc) = (0-5,10°, 4) (filled symbols) 
and (b) (^ 2 , 7 , J c ) = (0-5, 230°, 4) (empty symbols). The last column 
corresponds to the lowest energies for a given J after K— mixing 
j (Eq. 6) for the same intrinsic states as in (a) and (b). 


The decomposition and the energy E J K are quanti¬ 
ties that depend on the orientation of the principal axes of 
inertia with respect to the (x, y , z)-axes. Nevertheless, the 
following properties are deduced from the self-consistent 
symmetries imposed to the intrinsic states: = W^ K , 

E j k — E j _ k , and, if J is odd, the K = 0 component is 
forbidden. The dependence on K is removed once the K- 
mixing is performed since the norm W J and the energy 
Ep° j c are scalar quantities [10, 11, 19]. In addition, if 
J c = 0 the same probability distribution W J and the same 
angular momentum projected energy are found in the six 
sextants of the (/?2,7) plane. However, if the cranking term 
is non-zero, the (/?2,7) plane is split in two equivalent parts 
divided by the (7 = 120°, 300°) line. 

We now exploit these symmetries to perform consis¬ 
tency tests of the results and check the implementation of 
the method. Therefore, we select first two intrinsic states, 

| fa = 0.5,7 = 10°, J c = 4) and | fa = 0.5,7 = 230°, J c = 
4) that are symmetric with respect to the (7 = 120°, 300°) 
line (see yellow squares in Fig. 1(c)). We represent the de¬ 
composition of those states in components of the angular 
momentum J and intrinsic ^-projection K in Fig. 2(a)-(d), 
normalized to the total probability in a given J . Here we 
observe that the decomposition in K is different depending 
on the value of 7. For 7 = 10° the probability decreases 
in general rapidly with increasing K for a fixed value of J. 
Furthermore, the relative weight of the components with 
large K tends to increase with larger angular momentum 
J, while the K = 0 component for even J and K = 1,2 
components for odd J slightly decrease. These results are 
consistent with having the intrinsic long inertia axis nearly 
along the z- axis. On the other hand, the probability for 
a given J is distributed in a larger number of K compo¬ 
nents for 7 = 230° and these components are much flatter 
than in the previous case when the angular momentum J 
is increased. In this case, the intrinsic long inertia axis is 


almost oriented perpendicular to the z- axis. Nevertheless, 
the decomposition in J of both states, summing all of the 
K components, are identical. We observe two separate 
distributions for even J (Fig. 2(e)) and odd J (Fig. 2(f)), 
being the absolute scale larger for the former. The even 
(odd) distribution probability increases from J = 0 (1) 
until the maximum at J = 4 (5) is reached. Then, W J 
decreases, obtaining practically zero probability for even 
(odd) angular momenta larger than 16 (13). In Fig. 3 ) we 
represent the E^ifi 2,7, J c ) energies defined above. We see 
as in the previous case noticeable differences depending on 
the 7 values. For 7 =» 10° the energies rise rather quickly 
for large values of K while for 7 = 230° the energies are 
flatter. These differences are completely removed when K- 
mixing is performed through solving the HWG equations 
(Eq. 6) as it is shown in Fig. 3(c). There, three bands can 
be distinguished, namely, a ground state rotational band 
with A J = 2 built on top of J = 0^, and two A J = 1 
bands, being J = 2 + and J = 1 + the corresponding band- 
heads. 

We can test even further the performance of the angu¬ 
lar momentum projection by projecting the whole (^2,7) 
plane as it is plotted in Fig. l(d)-(f). For the J c = 0 case, 
the equivalence between the six sextants is preserved when 
angular momentum projection is performed. However, the 
angular momentum projected PES attained by restoring 
the rotational symmetry of the J c 7^ 0 states are symmet¬ 
ric only around the axis (7 = 120°, 300°) (Fig l(e)-(f)). In 
Fig l(d)-(f) only the PES for (J = J c ,cr = 1) are shown 
although the same equivalence is obtained for other values 
of (J, a). Apart from the symmetries discussed above, the 
angular momentum projection modifies significantly the 
surfaces obtained at the PN-VAP approach. In general, 
the minima found in the PES at the PN-VAP level are 
now wider and a slightly larger deformation is obtained 
whenever the angular momentum is restored. As a mat¬ 
ter of fact, these beyond mean-field correlations move the 
ground state from the spherical point to prolate configu¬ 
rations with /?2 ~ 0-5 (Fig. 1(d)), that was formerly a sec¬ 
ondary minimum in the PN-VAP calculation (Fig. 1(a)). 
This effect was already obtained with axial calculations [6] 
and is a self-consistent way to obtain the deformed ground 
state for the nucleus 32 Mg, i.e., as belonging to the ’island 
of inversion’ with an erosion of the N = 20 magic number. 

The last step in the present SC CM many-body method 
is the configuration mixing: 

\JM; NZ ; a) = ^ ffc„ tK>Je \JMK; NZ ; ft, 7) J c ) 

£ 2,7 ,K,Jc 

( 12 ) 

Again, the coefficients K Jc in Eq. 12, and the final 
spectrum, E Ja , are obtained by solving the general HWG 
equations: 

E (^W;K> - £ J WS'}) /£'} = 0 ( 13 ) 

{<*'} 

where {<7} = {^2,7, K , J c } now encodes all the constraints 
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Figure 4: (color online) (a) Excitation energies of the yrast states 
calculated for the nucleus 32 Mg with the GCM method with 17 ax¬ 
ial states and J c = 0 (Si), 49 axial+triaxial states and J c = 0 (S 2 ), 
81 axial+triaxial states and J c = 0, 2 and 113 axial+triaxial states 
and J c = 0,2,4. (b) Excitation energies for the three lowest states 
for each J-value, calculated with the GCM method with 113 ax¬ 
ial+triaxial states and J c = 0,2,4 (full symbols), large scale shell 
model calculations (open circles) and experimental data (asterisks). 


and K in a single index. T-L and M are the energy and norm 
overlaps respectively: 

•V [a}-{a'} = ( JMK ; NZ ’ ft. 7 , JcJMK 7 ; NZ- ft, 7', J') 

(14) 

n {«} Z { a '} = JcWMK'-NZ-^i, J') 

(15) 

To shed light on the impact of including time-reversal 
symmetry breaking states in the spectrum, the nucleus 
32 Mg has been computed with the GCM method using 
four sets of intrinsic wave functions. All of them are com¬ 
puted with nine major oscillator shells (7V 0 . S . = 9) in the 
working basis 2 . The simplest one (Si) is made of the 17 
axial and time-reversal symmetric states. Such states are 
marked in Fig. 1(a) with dots along the (7 = 0°,180°) 
axis. Then, the S 2 set is defined by adding 32 more time- 
reversal conserving states (J c = 0) in the (/? 2 , 7 ) plane 
(the remaining dots in the same figure). Finally, two more 
batches of states, S 3 and S 4 , are established by adding 32 
time-reversal symmetry breaking states with J c = 2, and 
32 more with J c = 4 (see the dots in Fig. 1). Therefore, 

S\ C S 2 C S 3 C S 4 , being the total number of states in 
the largest set equal to 113. 

The ground state bands calculated with the GCM method 
implemented with the different sets described above are 


2 We have checked the convergence of the energy spectra by per¬ 
forming axial calculations including up to thirteen major oscilla¬ 
tor shells. We have obtained practically the same results as the 
N 0 . s . = 13 ones already with N 0 . s . = 7, thus showing a good con¬ 
vergence of the results with respect to the size of the working basis. 


shown in Fig. 4(a). Firstly, the ground state energies ob¬ 
tained for the different calculations are pretty close except 
for the pure axial case, namely: 

E(0t) = -253.056, -253.477, -253.486, -253.498 MeV for 
Si, S 2 , S 3 and S 4 respectively. That shows that the ground 
state energy is converged with respect to adding time- 
reversal symmetry breaking components. However, the 
excited states are more affected by the inclusion of triaxial 
and cranking states. Hence, we first observe a moder¬ 
ate compression of the spectrum from the axial (K = 0) 
to triaxial calculations with J c = 0. The decrease in 
energy is larger with increasing the angular momentum, 
mainly due to the possibility of having more iT-mixing in 
the GCM states. However, the variational space for the 
excited states are much better explored if time-reversal 
symmetry breaking is allowed. Therefore, a significant 
compression of the spectrum is obtained for the S 3 and 
S 4 sets and the differences, once again, are bigger for 
larger angular momentum. In addition, we can infer that 
the excitation energies for the 2 + and 4+ states are al¬ 
ready converged with the S 3 calculation since they do 
not vary significantly from including J c = 4 states to the 
J c = 0,2 ones. This is not the case for larger values of the 
angular momentum, where probably intrinsic states with 
J c = 6 , 8 ,... should be also included in the GCM. For the 
sake of completeness, the full spectrum computed with the 
64 set is represented in Fig. 4(b). Here, the first two bands 
display a rotational character, with a parabolic trend in 
the excitation energies, 0 + band-heads and AJ = 2 spac¬ 
ing. A third band starting at 2% with AJ = 1 is also 
obtained, showing a slight odd-even J staggering. In ad¬ 
dition, large scale shell model (LSSM) results [24, 25] and 
experimental data [26-30] are also represented in Fig. 4(b). 
Thanks to the compression of the spectrum produced by 
the addition of cranking states, a remarkable agreement 
between the experimental and theoretical values for the 
2+ and 4+ energies is obtained. In addition, the present 
SCCM calculations predict very similar excitation ener¬ 
gies to the LSSM values for the g.s. band. However, the 
low excitation energy of the Oj state [30] is not repro¬ 
duced here. LSSM calculations have shown that this state 
is very sensitive to a subtle mixing of spherical Op-Oh and 
superdeformed 4p-4h configurations [25]. In the present 
framework, the inclusion of pairing fluctuations [9] and/or 
explicit quasiparticle excitations could help to solve this 
problem since the excited 0 + states are mainly affected by 
such a degree of freedom, lowering the excitation energies 
of those states. 

Finally, we explore systematically the effect of the in¬ 
clusion of time-reversal symmetry breaking states in the 
magnesium isotopic chain 24 - 34 Mg. The results are ob¬ 
tained with N Sm0 ' = 7 -the minimum that guarantees a 
good convergence with respect to the size of the basis in 
this isotopic chain- and the sets of wave functions defined 
above Si, S 2 , S 3 , he., axial and triaxial shapes with J c = 0 
and 2 are included. In Fig. 5 we plot the excitation en¬ 
ergies for the 2 + and 4+ calculated with these different 
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Figure 5: (color online) 2+ and 4/" excitation energies for the Mg iso¬ 
topic chain calculated with the GCM method including axial states 
(red squares), axial+triaxial with J c = 0 states (blue diamonds) and 
axial+triaxial with J c = 0,2 states (magenta open dots). Experi¬ 
mental values (black dots) are taken from Ref. [31] and references 
therein. 

approaches compared to the experimental values. We see 
that the axial calculations describe the trends of the ex¬ 
perimental data but the energies are largely overestimated. 
Including the triaxial degree of freedom without breaking 
the time-reversal symmetry reduces the excitation energies 
but the predicted values are still too high with respect to 
the experiments. Finally, adding J c = 2 states to the 
GCM set of wave functions compresses further the spec¬ 
trum and an outstanding agreement with the experimen¬ 
tal values is found. The only nucleus where the theoretical 
values tend to be lower than the experimental ones is the 
nucleus 24 Mg. Since this is a N = Z nucleus, some al¬ 
pha clustering and/or proton-neutron pairing correlations 
could be missing within the present framework which as¬ 
sumes a structure of the intrinsic states given by a direct 
product of protons and neutrons wave functions. However, 
mixing protons and neutrons to take into account such 
proton-neutron pairing correlations is beyond the scope of 
the present study. 

In any case, we have to underline that these results 
constitute the first explicit evidence of the compression 
of the spectrum when time-reversal symmetry breaking 
is taken into account in GCM calculations with particle 
number and angular momentum projection. Global cal¬ 
culations performed with these methods assuming axial 
symmetry have displayed a systematic overestimation of 
the excitation energies around a factor ^ 1.2 — 1.4 with 
respect to the experimental values, both for Skyrme [32] 
and Gogny [33] functionals. The present results show that 
such a disagreement could be corrected by including triax¬ 
ial and J c 7^ 0 states in the GCM framework. In fact, the 
incorporation of J c in the GCM ansatz (Eq. 12) is a gen¬ 
eralization of the double projection method of Peierls and 
Thouless [34, 35]. The double projection method is known 


to provide the exact translational mass in the case of trans¬ 
lations by taking as coordinates the position and the linear 
velocities in a generator coordinate method. We expect, 
therefore, that the moment of inertia of our theory will 
be similar to the one provided by the angular momentum 
projection before variation approach, instead of the Yoccoz 
moment of inertia given by the angular momentum pro¬ 
jection after variation method used in earlier approaches. 
This expectation is confirmed by our results that provide 
moments of inertia very close to the experimental ones. 

In summary, we have presented the first GCM calcu¬ 
lations with particle number and angular momentum pro¬ 
jection of HFB-like states considering different quadrupole 
deformations (axial and triaxial) and intrinsic cranking an¬ 
gular momentum. The performance of the method has 
been checked by taking advantage of the self-consistent 
symmetries imposed to the intrinsic many-body states. 
Since such wave functions were chosen to be eigenstates 
of a Dj h sub-group generated by the parity, simplex-x and 
T-simplex-y symmetry operators ({V, S x , <Sj}), the poten¬ 
tial energy surfaces (particle number and particle number 
plus angular momentum projected) must be symmetric in 
the (/?2,7) plane with respect to the (7 = 120°, 300°) axis. 
We have checked such a non-trivial property both in in¬ 
dividual states and in the whole (^2,7) plane, taking the 
nucleus 32 Mg as an example. The effect of including incre¬ 
mentally intrinsic states with more symmetries broken in 
the GCM framework has been also analyzed in 32 Mg and 
in the magnesium isotopic chain 24-34 Mg. The results 
have shown that adding (J c 7^ 0) time-reversal symme¬ 
try breaking states squeezes notably the spectra due to a 
better description of the excited states from a variational 
point of view. Such a compression puts the theoretical 
values on top of the experimental ones for the lowest 2 + 
and 4 + states in the chain. The next step will be the cal¬ 
culation of electromagnetic properties within the present 
SCCM approach and some work is in progress along these 
lines. 
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